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ABSTRACT 

A method for determining the free-field open-circuit voltage sensitivity 
of areversible underwater electroacoustic transducer from the difference 
in its input electrical impedance when loaded by water and air was 
investigated theoretically and experimentally. An equation for the 
Sensitivity was derived using reciprocal two-port network theory. The 
theory takes into account the diffraction due to the finite size of the 
transducer, its finite mechanical impedance, and its free-field radiation 
impedance in water. An experiment to test the predictions of the theory 
was performed using a 6" diameter, hollow, piezoelectric ceramic Spherical 
transducer. The result of the experiment agreed within several dB witha 
standard comparison calibration over the frequency range for which both the 
electrical impedance and comparison calibration data are considered 
reliable. The calibration method described, which has been termed the 
Delta-Z Method, may be useful for /a-s7tu monitoring of transducer 


sensitivity in installations which can be flooded and purged. 
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|, INTR it 


A. BACKGROUND 

For the purpose of the research described in this thesis, an acoustical 
tranducer is a device which transforms fluctuations of environmental 
pressure into electrical quantities or vice-versa. A reversible acoustical 
transducer is one which can act as either a transmitter or a receiver. 

There exist many calibration techniques which allow one to determine a 
transducer’s receiving sensitivity, that is, the electrical output for a given 
pressure field at the transducer face, and/or its transmitting response, that 
is, the pressure field resulting from a given electrical input. The focus of 
this study 1s on a method introduced by the Royal Australian Navy (RAN) 
whereby the open circuit voltage sensitivity of a reversible underwater 
transducer can be obtained from the difference in its input electrical 
impedance (AZ) when loaded by water and air. This method shall hereafter 
be referred to as the Delta-Z method. 

While there is agreement that the sensitivity can be obtained from AZ, 
there ts disagreement over the appropriate analysis of the Delta-Z method 
and consequently over the formula which applies (Ref. 1]. Preliminary 
investigations of the Delta-Z method applied to in-service transducers were 
"made by Taylor (Ref. 1] and Westbrook [Ref. 2]; however the experimental 
results were inconclusive. 

Because of its experimental simplicity, the Delta-Z method has great 


potential for n-s/tu transducer calibration. For this reason it was chosen 


for further study. It was hoped that, by using a simple geometry transducer, 
for which many of the needed parameters can be easily calculated, the 
controversy over the proper analysis of the method could be resolved and its 


usefulness established. 


B. OBJECTIVE 

The objective of the research described in this thesis is to develop, 
test and evaluate an underwater transducer calibration procedure for a 
simple geometry reversible transducer based on the difference in its input 
electrical impedance when loaded by water and air. 


C. QUTLINE OF REMAINDER OF THESIS 

In Chapter |! the various standard methods used for transducer 
calibration and the theory associated with the Delta-Z method are 
discussed. The experimental procedure and the raw data gathered from the 
spherical transducer are described in Chapter III, while the analysis of this 
data and its associated results are part of Chapter IV. A discussion of the 
results follows in Chapter V. A summary and some conclusions and 
recommendations appear in Chapter VI. 


I]. THEORY 


A. STANDARD CALIBRATION METHODS 

The calibration methods described in this section are discussed in 
greater detail by Bobber [Ref. 3] and some portions have been extracted from 
his book. | 

|. Comparison Calibration method 

In this procedure the electrical output voltage of the transducer 

under test is compared of to that of a standard transducer for which the 
sensitivity is Known. It is the most expedient and widely used method to 
determine the absolute sensitivity of most transducers. The results 
obtained using this method wil] be used throughout this thesis as the 
Standard to which the results of other methods will be compared. 

The open-circuit output voltage, @c, of the standard hydrophone is 
measured in a free-field position with its acoustic axis pointed towards the 
projector. The standard hydrophone is replaced by the unknown hydrophone 
and the open-circuit output voltage, ey, is measured. If the free-field 
voltage sensitivity of the standard is Mz, then the sensitivity of the 


unknown is found from: 


My = Mg ey . (Oa 
Cs 


2. Conventional Reciprocity Calibration Method 


A conventional reciprocity calibration requires three transducers 
of which one serves only as a projector (P), one is areciprocal transducer 
(T) and serves as both a projector and a hydrophone, and one serves only as a 
hydrophone (H). Any one of the three transducers can be the unknown or the 
one being calibrated; however, the calibration formula usually is derived for 
the free-field voltage sensitivity, My, of the hydrophone. The measurements 
are made in the far-field under free-field conditions. 

The arrangements and measurements made are shown schematically 


in Figure | and the procedure js as follows: 


a. the input current, ip to the projector is measured while the 


hydrophone free-field output voltage, epy, is also measured , 


b. the same process is repeated with the hydrophone replaced by the 


transducer and its output voltage, ep7, iS measured, 


c. the projector is then replaced by the transducer and both the 
input current, i7, and the output voltage of the hydrophone, eyy, 


are measured. 


INPUT HYDROPHONE OUTPUT 
CURRENT ae | VOLTAGE 


ie ——-] (P) =< (#) —a> "eu 


2) 'p — (°) — @ — & 


3 > —— (T) ————-) (4) —— °TH 


Figure 1. Diagram of the Measurements Needed for 
a Conventional Reciprocity Calibration. 


The free-field open-circuit voltage sensitivity, My, of the 


hydrophone is calculated using the following equation: 


My = | epyetH J | '/2 (2.2) 
epT'T 


where J is the spherical reciprocity parameter: 


J 22 de 
of 
In this case, d is the distance between the acoustic centers of the 


two transducers, p is the density of the medium and f is the frequency of 


operation. 


3. Iwo-transducer Reciprocity 
Only a single transmitting current and receiving voltage 
measurement (iz and ey, ) is required to perform a reciprocity calibration 
if two identical transducers are employed. In this case equation 2.2 is 


replaced by 
My = a i 1/2 (2.3) 


B. SELF-RECIPROCITY METHODS 

l. f —reci N 

A single transducer can be used in a two-transducer reciprocity 

calibration procedure if it is placed near a perfectly reflecting surface; 
reflecting a transmitted acoustic pulse back to the transducer, it receives 
its own transmission. This self-reciprocity arrangement is described by 
Bobber [Ref. 3] and is shown schematically in Figure 2 where the air-water 
interface is used as the perfectly reflecting surface. 


PERFECT REFLECTOR 
(WATER SURFACE) 


- 
oe 
“me 


<Ge Be 
| 
(x) TRANSDUCER 


Figure 2. Self-reciprocity Arrangement 


The resulting free-field open-circuit voltage sensitivity of the 
transducer, M.,, is determined by the expression: 


My = Ps - V/2 (2.4) 
IT 
wnere: 
eTy = open-circuit voltage across terminals during 
receipt of the reflected signal, | 
it = current flowing into the transducer during 


transmission, 


Js = spherical reciprocity parameter: 


Jo =20 (255) 
pf 


d = twice the distance from the radiating face of the 
transducer to the reflecting surface, 

p = the fluid mass density, 

f = frequency of the transmitted signal. 


2. pabin's Delta-2 Method 

Sabin (Ref. 4] pointed out that the ratio e¢,/i¢ was equal to the 
difference the input electrical impedance before and after receipt of the 
reflected pulse: 

AZ= 27 - Zp er, (2.6) 
IT 

where | 

Zy = input electrical impedance after receipt of the reflected pulse, 


Z-¢ = input electrical impedance before receipt of the reflected 


pulse, i.e. ina free-field. The expression for the open-circuit 


sensitivity from equation 2.4 then becomes 


My = (AZe)* (2..7) 


3. Royal Australian Navy Method 


The Royal Australian Navy (RAN) developed a calibration method, 
claimed to be an extension of Sabin's method, in which the depth of water 
for the reflected impedance measurement is taken to be zero (the 
transducer is in air) [Ref. 5]. The planar reciprocity parameter is used to 
determine the sensitivity instead of the spherical reciprocity parameter. 

The electrical impedance of the transducer to be calibrated is 


measured both in water (Zwe¢) and in air (Zqry) where: 
AZ = Zwet a Zdry ’ (2.8) 


so that the open-circuit sensitivity now becomes: 


My = (AZ Jp) 1/2 (2.9) 
where: 
Jp = the planar reciprocity parameter, . 
Jp =2A , (2.10) 
pc 


A = area of the active face of the transducer, 
p = density of the medium, 
C = speed of sound propagation in the medium, and 
My = free-field voltage sensitivity. 
The method described above will be hereafter referred to as the 
Australian or RAN method, while the method that will be developed in the 
following section will be referred to as the Delta-Z method. 


C. RECIPROCAL TWO-PORT NETWORK ANALYSIS OF THE DELTA-Z METHOD 

The treatment of the Delta-Z method as an extension of Sabin’s method 
is not justified. For one thing, in conventional reciprocity methods each 
transducer is assumed to be in the far field of the other. Applied to 
reflection methods, this means that the transducer must be in the far field 
of its image in the reflecting surface. This condition cannot be met for a 
finite size transducer as the distance to that surface approaches zero. 
Furthermore, as will be apparent from the theory developed in this section 
and in the experimental results presented in Chapter V, the appearance of 
the planar reciprocity parameter in equation 2.9 is incorrect 
(notwithstanding that equation 2.9 does not include the effects of 
diffraction and finite mechanical impedance). 

The proper analysis of the Delta-Z method is in terms of a reciprocal 
Cwo-port network. Such an analysis is given in this section and leads to a 
different equation for the open-circuit sensitivity than equation 2.9. 

A reversible electroacoustic transducer may be represented as a 


reciprocal two-port network, as schown in figure 3: 


Electrical Transducer Mechanical 


Port Port 





Figure 3. Two-Port Representation of an 
Electro Acoustic Transducer 


The appropriate two-port equations are 


Bare e ie ew, (Zip 
p=Ti+2yU, (23 
here: 
e = the voltage across the electrical terminals, 
| = the current passing through the electrical terminals, 
U = the volume velocity directed into the face of the transducer, 
p = the average acoustic pressure at the face of the transducer, 
Zp = the electrical impedance with blocked mechanical terminals, 
Zr = the mechanical impedance with blocked (open-circuit) electrical 


terminals,and 
T= the transduction coefficient for a reciprocal transducer. 
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Now consider such a transducer radiating into air as schematized in 
Figure 4. 


=_- 


Electrical Transducer Mechanical 
Port Port 


p=0 





Figure 4. Transducer Radiating in Air 


The radiation impedance of the transducer in air is considered to be 
negligible compared to Z,, and hence is represented in Figure 4 by a short 


circuit, 
p=0. (2.13) 


Solving for U in equation 2.12 with p=0 and substituting in 2.11: 


eal, 
Lm 


eeeeer | Sala ti) 
Lin 


Thus the electrical input impedance in air becomes: 
Zaip=ze - 12 @ (2.14) 
Lm 


Now consider the transducer transmitting into a free-field in water, as 


represented by Figure 5. 


Electrical Transducer Mechanical Acoustic 
Port Port Field 





Figure 5. Transducer Radiating in Water 
Here Zrag rp is the free-field radiation impedance of the transducer in 
water. Hence: 
DP = ~Zeag pe U . (2.15) 


Substituting this into equation 2.12 gives 


-Leag fp U= Tit 2, U 
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or 
Us __-T | 
(Zm + Zrag ff) 


substituting this into equation 2.11 yields 


e=Zci+T Ti } 
Lm + Zrad ff 


The input impedance therefore becomes: 


2 
aie ee- [7 (2.16) 
Lm * Zrad ff 


and the impedance difference between water and air is given by: 


Al=Zwater-Zair=Ze- — 1%  - Zee T? 
Lm * Zrad ff Zin 


which simplifies to 


AZ= oT *Zrag ee (2.17) 
(Zr + Zrag tt lm 


Thus the electromechanical coupling coefficient T is proportional to the 


Square root of the difference in the input electrical impedance of the 


transducer when loaded by water and air. 


ZI 


Note that if Z-aq ¢¢ and Zp, are Known, then T may be found from AZ 
using equation 2.17. If Zp is also Known, then all the two-port parameters 


would be Known and in principle one may analyze any proposed transducer 


application. In many cases of interest, however, Zr, and Zp need not be 
known to obtain the sensitivity of the transducer, Mo, provided that [azI/I2! 


is significant. As will be shown below, this simplification is obtained when 
the radiation impedance is smal] compared to the mechanical impedance, 
which is very often the case for piezoelectric ceramic transducers except 
near frequencies of mechanical resonance. This is the great utility of the 


Delta-Z method. 
To relate the free field open circuit voltage sensitivity, Mo, to AZ, 


imagine an acoustic wave impinging on the transducer. The network 


representing this situation is depicted in Figure 6. 


Mechenical Acoustic 
Field 


Electrical 
Port 


Trensducer 


Port 
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Figure 6. Transducer in an Acoustic Field 
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Where: 
Dp = the blocked pressure, that is, the average pressure 
over the face of the transducer with zero 
volume velocity, 

Dre = the free-field pressure, that is, the pressure which 
would exist at the location of the transducer if it 
were absent, 

D = the diffraction constant, and 
Zeag tf = the free-field radiation impedance of the 
transducer in water. 


By definition, the blocked pressure and the free-field pressure are 
related by the diffraction constant, D, which depends only on the geometry 
of the transducer and on the acoustic wavelength: 


Db =D Der. (2.18) 


The free-field open circuit (i=0) voltage sensitivity is derived from 
Figure 6 as follows. 


The free-field open-circuit voltage sensitivity may be cast as 


Moc ff" Mo™ foc “foc P Po: “ 
Prf DP Ph Prt 


When i=0 is substituted in equations 2.11 and 2.12, they become: 


aS 


p= U Lm 
giving: 
/ Coc = Ts. (2.20) 
p Lm 

Also, from Figure 6, 

D= Py- Zragrf UY (2.21) 
Therefore: 

pil + Zrad tt] = Pb = 2m * Zrad ff 

Zm Zm 

giving: 


D= cm | (2.22) 
Pp = 2m * Zraad rf 


Substituting equations 2.18, 2.20 and 2.22 into 2.19, the open circuit free- 
field sensitivity becomes: 


My = 7. | Zm | D (2.23) 
f 


Lm |2m* Zrad f 


The transduction coefficient T is found from equation 2.1 7: 


T =/AZ Lin (Zr) u Lead ff) W/2 
crad ff 


Equation 2.23 then becomes: 


Mo =|42 Zm (2m * Zrad tr|'”2 D Zin 
Zrad ff Lm (Lm + Zraa ff? 
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which simplifies to 


My =D ee V2 (2.24) 
Zrad ff ‘2m * Zrad ff? 


which is the desired result. 


If Zrag tr << 2mm, aS IS USUally the case for a piezoelectric transducer 


except near frequencies of mechanical resonance, equation 2.24 becomes 


much simpler: 
My =D | _az V2 (2.25) 
Zrad ff 


Note that in this case the mechanical impedance need not be known to 
determine Mp. This feature makes the Delta-Z method potentially very 
useful for n-s/tu calibration of in-service transducers. Also note the 
Similarity and difference between this and the Australian result for the 


Sensitivity which was given in section I1.B.3 and is reproduced here: 
My = (AZ Jp) UZ 


It is important to note that Zraq ff in equation 2.25 is not a reciprocity 
parameter and should not be construed as the inverse of the planar 
reciprocity parameter, Jp, even if equation 2.9 took diffraction into 
account, which it does not. These differences will be apparent in Chapter V 
where the experimental results of both methods will be compared to the 
standard comparison calibration of a spherical transducer performed at 
TRANSDEC, San Diego. 
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A similar two-port analysis may be performed which retates the short- 
circuit free-field current sensitivity of a reversible transducer fo the 
difference in the input electrical admittance when the transducer Is loaded 


by air and water. The result is 


Moc = D} AY Yraare | '. (2.26) 
(1+YmZrad ff) 


where: 
Moc = the short-circuit free-field sensitivity, 
AY = the change of admittance between the water and air, 
Yrad ff = the radiation admittance of the transducer, 
Yo = the mechanical admittance with shorted electrical 
terminals, and 
Ye = the electrical admittance with shorted mechanical 
terminals. 
This formula may be more appropriate for a low mechanical impedance 
device such as a moving coil transducer. (Note that a moving coil transducer 
is anti-reciprocal and that we must replace T@ in the above by -T2. This has 
no effect on the magnitude of the final result.) . 


In the case where 


Ym 2rad ff << |. 


then equation 2.26 simplifies to 


26 


The open-circuit voltage sensitivity can be determined using 
admittance parameters instead of impedance parameters with the following 


result: 


Ye + (YE Ymn~ T2)Zpad tf 


where: 


T = the transduction coefficient for a 


reciprocal transducer, and 


T = LAY (Ym*¥rag ep) 72 
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Il. DESCRIPTION OF EXPERIMENT AND RAW DATA | 


A. DESCRIPTION OF THE TRANSDUCER 

The transducer used for this experiment was a 6-inch outer diameter, 
hollow, piezoelectric sphere. It was borrowed from TRANSDEC and was 
built at NOSC, San Diego. It consists of two radially poled piezoelectric 
ceramic hemi-spheres of unknown (probably about 3/8") thickness glued 
together. It is potted in polyurethane of approximately 1/4 inch thickness 
and has a 12 inch long hard rubber handle to house the connection to the 100 
foot long electrical cable (Cable * 71102, DSS-3 (RF), MIL-C-915A), as 
depicted in Figure 7. Its low frequency capacitance is 210 nF (@1000 Hz). 


B. EXPERIMENTAL SET-UP AND PROCEDURES 

Measurements of Zwater and Zajp were carried out at TRANSDEC, NOSC, 
san Diego. The main feature of this facility is a large water tank which 
Simulates free-field conditions for a transducer placed at a certain position 
(6m depth at the center of the pool) (Ref. 6]. Its unique construction focuses 
the sound energy on the outside perimeter of the tank where it is reflected 
into a sound trap, eliminating the energy reflected back to the transducer. 

Measurements were conducted using a continuous wave (CW) signal 
generated from the HP-4192 Impedance Analyzer which was controlled by a 
HP-9826 computer and HP Basic 4.0 software. The equipment set-up of 
Figure 7 was used to determine the resistance (R) and reactance (X) 
components of the impedance (Z) of the spherical transducer in both water 


and air. 
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HP 41924 Low Frequency 
Impedance Analyzer 





Figure 7? : Equipment Configuration Used to Determine 
sensitivity by the Delta-Z Method . 
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The procedure was as follows: 

The sphere was placed 6 m below the surface of the water for half an 
hour to allow its temperature to reach equilibrium with the surrounding 
water. The Rwater and Xwater Components were measured and stored on a 
floppy disk. The data covered the range of 1.0 KHZ to 15.0 KHz in 100 equal 
increments. The sphere was then pulled out of the water, and Raj, and Xair 
were immediately measured at the same frequencies as in water and stored 


on the same file as the water data. 


C. RAWDATA 
Figures 8 to 11 show the raw impedance data obtained for the sphere 


following the above procedure. 
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R water (Ohms) 


X water (Ohms) 


R water vs Frequency 
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Frequency (kHz) 


Figure 8. Rwater aS a Function of Frequency. 
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Figure 9. Xwater as a Function of Frequency. 
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X air (Ohms) 


R air (Ohms ) 
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R air vs Frequency 
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Figure 10. Rajp as a Function of Frequency. 


X air vs Frequency 
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Figure 11. Xajp-as a Function of Frequency. 
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IV. ANALYSIS AND RESULTS 


A. EXTRACTION OF AZ FROM RAW DATA 
The electrical impedance measured at the terminals of the transducer 


is represented by the complex quantity: 


Z=R+jX 
where : 
Z = impedance, 
R =resistance, and 
X = reactance. 


For the measurements taken in water the expression is written: 
water = Rwater * J Xwater 
and similarily for the air measurements: 
Zair = Rair * J Xair. 


The magnitude of AZ for insertion in the sensitivity calculations is given 


Dy: 
1/2 
| AZ| = neers. as Rajr ) 4 + (Xwater a Xair ) 2] (4.1) 


and the angle associated with it is: 
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{Delta-Z'(/Z water 


ZAZ=arc tan Xwater - Xair (472) 


Rwater ~ Rair 


Only the magnitude of AZ was used in this experiment. 
The values of resistance and reactance used in equation 4.1 were 
obtained as described in section 111.B. Figure 12 shows a plot of 


IAZI/IZwaterl VS frequency. 


[Delta-Z!/Z water vs Frequency 
.42 
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Figure 12. Plot of IAZI/IZwaterl vs Frequency 
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B. VARIABILITY OF AZ MEASUREMENTS 


To ascertain the reliability of the values of AZ, the experiment was 
repeated three different times. Figure 13 depicts the results obtained for 


IAZI/IZywater! from these three independent experiments. 


LDelta-Z)J/“7Z water vs Frequency 
-oo 


a 
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Figure 13. Plot of IAZI/lZwaterl for Three Independent 
Measurements 


The results obtained indicate that random experimental error is 
probably less than 1% for frequencies below resonance. A detailed error 
analysis is included in Appendix A. Determination of the remaining 


parameters needed to calculate Mo (D, Zpaq rr, 2m) is described below. 
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C. DETERMINATION OF D, Zrag rp and 2m 
|. Diffraction constant (D) 
Recall from Figure 6 and the discussion which followed that the 
diffraction constant is defined to be the ratio of the blocked pressure to the 


free-field pressure. It was derived by Henriquez for a sphere [Ref. 7] and is 


given by : 
Dee ree (4.3) 
\/1 + (kr) 2 
where: 


D = magnitude of the diffraction constant, 
kK = wave number = a/C, 


a = radius of the sphere. 


The diffraction constant of a transducer with a more complex 
geometry is not necessarily a simple expression and may have to be found 
experimentally. Figure 14 depicts the diffraction constant as a function of 


frequency for several interesting shapes. 
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Figure 14. Diffraction Constant as a Function of 


Frequency [Ref. 7] 


2. Radiation Impedance (Zraq rf) 


The radiation impedance of a sphere in radial motion is given by 


Beranek [Ref. 8]: 


Z rad ff = Radiation impedance In free-field 


= — Japgca 
S (c + joal! -jcd)) 
@) 


(4,4) 


where: 
w = radial frequency, 


Pg = density of the medium, 
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c= speed of sound propagation in the medium, 
S = surface area of the sphere, 
= radius of the spherical transducer, and 
= absorbtion constant. 
Since c/w is in the order of 1x1079 for all frequencies of interest in this 
study, it can be neglected without significantly affecting the results [Ref. 
8, pp 57,58]. Therefore it is assumed: 


s [1+(ka) 2] 


The real and imaginary parts of equation 4.5 are : 


Rrad = Poc(ka)* — = (ka) Xpag (4.6) 
S (1+(ka) 2) 
S (1+(kr) 2 
3. Blocked Mechanical Impedance (Z,,) 


The blocked mechanical impedance of the sphere was extracted 
from its input electrical impedance in air and water as described below. 
The success of the technique employed depends upon the assumption that the 
radiation resistance in water is significant compared to the blocked 
mechanical resistance. 

The input electrical impedance of a simple reversible piezoelectric 
ceramic transducer in air may be modeled by the approximate equivalent 


circuit shown in Figure 15 (Ref. 9]. 
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Figure 15. Approximate Equivalent Circuit of a Simple 
Piezoelectric Transducer in Air. 


where: 


Cg = open circuit voltage, 

Cy = blocked electrical capacitance, 

Cm = blocked mechanical capacitance, 
Ly = blocked mechanical inductance, and 
Rm = blocked mechanical resistance. 


The notation follows Wilson [Ref. 9], i.e. capital letter subscripts refer to 
equivalent electrical quantities, lower case subscripts refer to mechanical 
quantities. Only one mechanical resonance has been modeled in Figure 15; in 
principle any number of additional resonances can be modeled by adding one 
more paralle] R-L-C branch for each. 

In making the above electrical analogy it is assumed the vibrating 
elements are mechanically equivalent to a single degree-of-freedom mass- 
spring-dashpot system and that the radiation impedance in air can be 


neglected. It is also assumed that the equivalent circuit parameters are not 
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frequency dependent. The latter assumption is not crucial to the success of 
the Delta-Z method, however, since the value of the mechanical impedance 
significantly affects the calculated sensitivity only near frequencies of 
mechanical resonance. 

When the transducer is surrounded by water the radiation impedance 


cannot be neglected and is represented in Figure 16 by Rpap and Lpap. 


Cy 
8p Ru * Rrap 


Ly + LRap 


Figure 16. Approximate Equivalent Circuit of a Simple 


Piezoelectric Transducer in Water. 


The procedure used to extract the mechanical impedance of the 


sphere, or, more precisely, the mechanical resistance Rp, inertance Lppy, and 
compliance Cr, is as follows. 

At the frequency of mechanical resonance in air (wp), the effects of 
Lm and Cy cancel and leave Ry in parallel with Cy (see Figure 15), sie the 


following expression for the input electrical impedance: 


I +(WpRuCo) 2 


where : 


@,= the natural mechanical resonance frequency. 


Separating Zajr into its real and imaginary parts we obtain: 


Rair = Ry 
1 +(@gRmCo)? 
1+(WpRvCo)2 


Rearranging equations 4.9 and 4.10 gives: 


Nair =~ Gp Ry Co. 
Rair 
Rair® *Xair® = Re * 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


Values of Raip and Xa; values are obtained from the input electrical 


impedance data with the transducer in air. Rr and Cp are then easily found 
by inverting equations 4.11 and 4.12. The same procedure is used to find the 


values of (Ry + Rrap) and Cy with the transducer in water (Figure 16). The 


values Of Ru, Co, and Raap at resonance found this way are: 
Ru = 1290 
Cy = 209.nF 
RR Ap =1.89x 1079. 
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Using the fact that the transformation of the mechanical 
resistances, Rr and Rrag, into their equivalent electrical quantities, Ry and 
Reap, 's identical, the following may be written: 


Ry = RM Reaq - (4.13) 
RRAD 


Hence Rp, may be found from this equation using the values of Raap and Ry 
extracted previously once Rrag iS known. Since Rag can be easily 
calculated for the sphere, this is a simple procedure. The determination of 
Read Will be one of the crucial factors on which the successful application 
of the Delta-Z method to 7-s/tu transducer calibration will depend. (It is 
not necessary to know Rrag, however to monitor degradation in sensitivity.) 

It should be noted that the ratio Ry,/Rrag is related to the 
electroacoustic efficiency of a piezoelectric transducer, n,as defined by 
Wilson (Ref. 9,p. 39] and Hunt (Ref. 10], by: 


Ry = Rm =lon (4,14) 
RRAD Rrag 0 
So that equation 4.13 may be written 


Rm = 1-1 Rrag 
n 
Having found R~, from equation 4.13, Lp, and Cy, are found from the 


center frequency, @), and quality factor, Q, of the mechanical resonance 
using 


Lm = QRm (4.15) 
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and 


Oy QRm 


The resulting values of Rm, Ly and Cp, for the sphere are: 


Rm = 5.93 x 106 Nsec/m®? 


Lin = 1.25 x 103 Nsec/m® 


Cry = 1.46 x 10-15 m/N sec, 


and the electrical values of the electrical equivalent values Rwy, Ly and Cy 


are Ry = 129 0 
Lm = 27.2 mH 
Cy = 6.7 nF, 


Figure 17 shows a plot of Z,, calculated using these values to the 
radiation impedance calculated for the sphere. 
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Figure 17. Zp and Zpag pr as a Function of Frequency 


V. DISCUSSION 


A. COMPARISON OF DELTA-Z METHOD TO STANDARD CALIBRATION METHOD 

The standard of comparison for all other calibration methods used to 
determine the sensitivity of the spherical transducer is the comparison 
calibration described in Section II.A.1. Tne comparison calibration of the 
spherical transducer was obtained at the Transducer Evaluation Center 
(TRANSDEC) at the Naval Ocean Systems Center (NOSC), San Diego, Ca, using 
a Navy standard type F-37 transducer. The comparison calibration method is 
widely used by the Navy and most other organizations that require absolute 
calibration of hydrophones. 

Figure 18 shows the results obtained at TRANSDEC by the Delta-Z 
method using eq. 2.24 (curve labeled “Delta-Z Method") compared to those 
obtained from the standard comparison calibration method. The two 
methods give results which agree within several dB except at the extremes 
of the frequency range of the experiment. The results from the Delta-Z 
method should be disregarded for frequencies above approximately 12 kHz 
Since the second and third mechanical resonances occur in this region and 
were not accounted for in the estimation of Zp). Also, the results from the 
comparison <alibration are suspect at frequencies below several kHz, as the 
pool at TRANSDEC cannot be considered anechoic at these low frequencies. 
In the frequency range 3 to 12 kHz the largest discrepancy between the 
results of the two methods is approximately 3 dB while the results over 


most of the frequency range agree within about | dB. 
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Figure 18. Sensitivity Obtained From Various Calibration Methods 


The importance of including the effects of finite mechanical impedance 
in the Delta-Z method in order to obtain accurate results near frequencies 
of mechanical resonance is readily apparent. The curve labeled “Without Z,," 
in Figure 18 was calculated using equation 2.25, which assumes 
Zm?>Zrad ff. The error made at 11 kHz by not accounting for finite Z, 


amounts to more than 5 aB. 


The close agreement between the results of the Delta-Z and comparison 
calibration methods below the frequencies of mechanical resonance was 
quite unexpected and js very encouraging for the prospect of applying this 
method to the monitoring of in-service transducers. Before the Delta-Z 
experiment was performed it was anticipated that the precision of the 
results would be best for frequencies around mechanical resonance, as Zp is 
minimum there and the effect of water loading would therefore be more 
observable, i.e. IAZI/|Z| would be larger. It was for this reason that a 
relatively large (by laboratory standards: 6" diameter) sphere was chosen 
for the experiment. Based on the experimental results there is good reason 
to believe that if the diffraction constant and the free-field radiation 
impedance can be determined within a few percent, the sensitivity of an in- 
Service transducer can be determined within | dB if AZ can be determined 
within 2 dB. From the magnitude of AZ/Z observed in this experiment 
(typically a few percent, see Figure 12), the precision with which 


impedance needs to be measured is probably about 0.1%. 


B. COMPARISON OF THE ROYAL AUSTRALIAN NAVY METHOD TO STANDARD 
CALIBRATION METHOD 
As was described in section I1.B.3, the equation resulting from the 


Royal Australian Navy analysis of the Delta-Z method is: 
My =(A/ Jp yi/2 


where: 


My = open-circuit sensitivity of the transducer, 
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AZ = the change of input electrical impedance between 
water and air, and 
Jp = the planar reciprocity parameter. 


The results obtained using this equation are labeled “Australian Method” 
in Figure 18. The sensitivity calculated using the RAN method is in great 


disagreement with. that calculated using the other two methods. 


MMARY, CON ION AND R MENDATIONS 


A. INTRODUCTION 

The objectives of the research described in this thesis were 1) to 
develop an underwater transducer calibration procedure based on the 
difference in input electrical impedance when loaded by water and air, and 
2) to experimentally test the procedure using a simple geometry (spherical) 
transducer. 


B. SUMMARY 

Using reciprocal two-port network theory, an equation for the free-field 
open-circuit voltage sensitivity of areversible underwater acoustic 
transducer was derived in terms of the difference in input electrical 
impedance when loaded by water and air. The theory takes into account the 
diffraction due to the finite size of the transducer, its finite mechanical 
impedance, and its free-field radiation impedance in water. An experiment 
to test the predictions of the theory was performed using a 6° diameter, 
hollow, piezoelectric ceramic spherical transducer. The results of the 
experiment agreed within several dB with a standard comparison calibration 
over the frequency range for which both the Delta-Z analysis and the 


comparison calibration data are considered reliable. 
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C. CONCLUSIONS 

It has been established that the open-circuit free-field voltage 
sensitivity of a reversible underwater piezoelectric acoustic transducer can 
be obtained by the so-called Delta-Z method, that is, from the difference in 
its input electrical impedance when loaded by water and air. The appropriate 


formula, derived from reciprocal two-port network analysis, is given by 


Mo =D AZ Zp 1/2 
Zrad ff “4m * Zraq ff? 


The finite mechanical impedance of such a transducer must be taken into 
account in order to obtain precise results near frequencies of mechanical 
resonance. Apart from such frequency regions the radiation impedance is 
generally negligible compared to the mechanical impedance and the 


sensitivity can be found from the approximate formula: 


Mo=D[ az | ¥2 
Zrad ff 


In this case neither the blocked impedance nor the open-circuit mechanical 
_ Impedance need to be known to determine the open-circuit sensitivity. 

The Delta-Z method should not be regarded as a high precision method -- 
the precision of the results rests on the precision with which the difference 
in two quantities which often are nearly equal can be measured. 
Nevertheless, because it is such a simple procedure, it may be useful for in- 
Situ monitoring of in-service transducers in installations which can be 


flooded and purged. 
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| RECOMMENDATIONS 
Recommendations for further study are: 

1. Investigate the Delta-Z calibration method using piezoelectric ceramic 
transducers of more complex geometry with the aim of applying this 
method to the 7a-s/tu monitoring of in-service transducers. 

2. Investigate the application of the Delta-Z method to a low mechanical 


impedance transducer such as moving coil transducer 
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APPENDIX 


ERROR ANALYSIS 


A. GENERAL 

Three types of errors are considered in this analysis: 
|. Random error, which is an indication of the repeatability of a 

measurement. This was determined experimentally by repeating the AZ 
measurement three times as described in Section IV.B. 

2. Estimated calibration error, which is an approximation of the 
systematic error of each measurement. 

3. Approximation error, values for the different parameters used in the 


calculation of sensitivity had to be truncated and/or estimated to a 
certain precision. 


B. ERROR IN AZ 


1/2 
[AZ| = | (Rwater-Rair!? (Xwater-Xair!”. 
where: 


AZ = change of impedance for the transducer from water into air, 
Rwater = Resistance of the transducer in water, 

Rajr = Resistance of the transducer in air, 

Xwater = Reactance of the transducer in water, 


Xair = Reactance of the transducer in air. 
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Using standard statistical analysis, the standard deviation for the Az 
parameter is: 


SA7 = 2 2 2 2 
= : a ; SRwater! ( =) SRaip * 
AZ dRwater SRair 


1/2 
2 2 2 2 

: 1 ‘SXwater’ “* —| ‘Sain? 

d&Xwate 8Xair 


The above assumes that Rwater, Rair, Xwater and Xair are statistically 








independent. 

lf we further assume that 

Rater . Ro ce > oaatan = Xai = 5, then the fractional deviation of 
AZ is given by 


1/2 
Saz_| 25? Rwater ~ Rair? * 2S? water ~ Xair” 7 


AZ (Rwater ~Rair)? *  (Xwater ~ Xair! 


which leads to 


oe-* 
AZ AZ 
From the HP-4192A handbook; S$ = 0.3% of reading, so that 
Saz = 0.6% 
AZ 
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C. ERROR IN THE BLOCKED MECHANICAL IMPEDANCE, Zr, 
Lm = Rm u jl@Lrpy v 1/@Cy,) = Ry + | Xm 
where w is the radial frequency and L,, and C,, are given in terms of Rp», by 


Lr, = QRp and I/Cmy = Wy QAR, . 
Wo 
where: 


Wo = mechanical resonance frequency (rad/sec) 


Q = quality factor of mechanical resonance. 


Then 


Slim =Q8Rm = and = 8(1/C§p) = QD ERey 
Wo 


ignoring the errors in f,@) and Q, compared to the error in Rp, 
then, 


SZ = 1 +j( 2nfQ+ wQ) bRry, 


W 2nf 
and 
tm ~ Rm 
fm = Rm 
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Assuming an estimated 10% error on the blocked mechanical resistance Ry, 
(from the estimated 10% error on the equivalent electrical resistance Ry), 


ey 10% 


Zm 


It should be noted that the influence of errors in Zp, on the error in Mo 
will only be significant near resonance, where Zp, is of the same order of 


magnitude as Zrag ff. 


D. ERROR IN THE RADIATION IMPEDANCE, Zpag rf 


For the spherical transducer used in this project: 


S [1 + (ka)2] 
IZpaq ffl = Poe Ka 5 (pC) | + (ka)~2 | a 
S [1 + (ka)2] 1/2 S 


taking the log of both sides and differentiating: 


SlZrag ffl = SPC - 85 + Ska) _ | 


rad ff! poe 9S (kadS. [1+ (ka)? ] 


SEpDe poo * oka) | 
po¢ «6S tka)~—sd 1 + tka? ] 


5 


but S28 = -2 64 = -2 S(ka) ignoring the error in k 
ahd a 
[1 + (ka)?] 


so that | -O2 * -&ka) |] 
Ss (ka) [1+ (ka)?] 


lies between Ka) and 2 aka) 


Ka Ka 
then 
? 2| 1/2 
“Zrad tt € |Spoc * 2 Skad 
cradff | Pos Ke 


= 2 S(ka), ignoring the error in PoC. 
Ka 


The fractional error ina, therefore in ka, is estimated to be about 5%. 
hence é 


“IZpag ffl ~ 10% 
IZead ff! 
E. ERROR IN THE DIFFRACTION CONSTANT, D 


Di | ! 


Vi + (ka)2 
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D 1+(ka)? ka 
F, ERROR IN THE SENSITIVITY, Mo 


From Chapter II, 


7 1/2 
Mo =O 2» | (Qoay 
Zrad tf ‘Zrad tf * 2m) 


For simplicity, define []= AZ Lm 
Zrad ff “Zrad tt * 2m) 


and Zrag f= Zr. 


Differenciating equation 2.24, 


dMo = My 2 * ANDI (alee! Xx 
Dy 2 


AZ Zp? tZnZr ( Z 4 +Zp Zp)? 


which can be written as 


8Mo = = §D +D AZ + Zr (62/2) + (22, +Zry M82p/ 2p) 
Mo DB 2 AZ (Z,. +Zm) 


but from sections D and E above, 


8D <&(ka) and 82, ¢ 2 &ka) 
D <) Zr Ka 


So that 


a7 


8M = 1+D(2Z-+2m) (ka) * D BAZ+ Zp 82m 
Zp +2Zm) ka 2 AZ (Zp+Zm) Um 


Near resonance take 2. =Z,, andD< 1. 
Then 

SMo < - . 7 CC ie 
Mo a 


Smo = [(2 x 5%) 2 + (0.5 x 0.4%) 2+ (0.5 x 10%) 2]? 
Mo 





SMoz 11% 
Mo 


Away from resonance let Zp << Zp, and let (82m/2m) ~ (8Z-/2-) and take 
D<1, so that 


ao = 8D + D Baz + 82p 
Mo D 2 A 2, 
= 28&ka)+ 1 8aZ 
ka 2 AL 
SO 
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SM =| /2 Stka)\24/ Saz\ *| 1” 
0 ka 2 AZ 


SMy< [ (2x52)? + (05x 0.4%)2] |” 


SMo < 10% . 


Converting the result to decibels: 


Error in dB = 20 log | +( SMo/Mo? = 20 log 1.1 = 0.8 cB 


Thus the estimated calibration error is approximately */. 1 dB for all 
frequencies. 
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